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ENVIRONMENTAL DEPENDENCE OF GALAXY MERGER RATE IN ACDM UNIVERSE 

HUNG-YU JlAN 1 , LlHWAI LlN 2 , AND TZIHONG ChIUEH 1,3,4 

ABSTRACT 

We make use of four galaxy catalogs based on four different semi-analytical models (SAMs) imple- 
mented in the Millennium simulation to study the environmental effects and the model dependence 
of galaxy merger rate. We begin the analyses by finding that galaxy merger rate in the SAMs has 
mild redshift evolution with luminosity selected samples in the evolution-corrected B band magnitude 
range,— 21 < Mg < —19, consistent with results of previous works. To study the environmental de- 
pendence of galaxy merger rate, we adopt two estimators, the local overdensity (1 + S n ) defined as the 
surface density from the n'^-nearest-neighbor (n = 6 is chosen in this study) and the host halo mass 
Mh- We find that galaxy merger rate F mg shows strong dependence on the local overdensity (1 + S n ) 
and the dependence is similar at all redshifts. For the overdensity estimator the merger rate F mg is 
found about twenty times larger in the densest regions than in under-dense ones in two of the four 
SAMs while it is roughly four times higher in the other two. In other words, the discrepancies of the 
merger rate difference between two extremes can differ by a factor of ~ five depending on the SAMs 
adopted. On the other hand for the halo mass estimator, F mg does not monotonically increase with 
the host halo mass Mh, but peaks in the Mh range between 10 12 and 10 13 M , which corresponds 
to group environments. High merger rate in high local density regions corresponds primarily to the 
high merge rate in group environments. In addition, we also study the merger probability of "close 
pairs" identified using projected separation and line-of-sight velocity difference C mg and the merger 
timescale T mg ; these are two important quantities for observations to convert the pair fraction N c 
into galaxy merger rate. We discover that T mg has weak dependence on environment and different 
SAMS, and is about 2 Gyr at z ~ 1. In contrast, C mg depends both on environment (declining with 
density) and different SAMs; its environmental dependence is primarily due to the projection effect. 
At z ~ 1, it is found that only ~ 31% of projected close pairs will eventually merge by z — 0. We 
find that the projection effect is the dominant factor to account for the low merger probability of the 
projected close pairs. 

Subject headings: galaxies: evolution — galaxies: interactions — cosmology: theory 



1. INTRODUCTION 

In a standard A cold dark matter (ACDM) model, 
the major framework is the hierarchical structure for- 
mation that dark matter halos and galaxies within 
them are assem bled from the successive accret i on of 
smaller objects ([Peebles! 119821: iBlumenthal et all 11984 
iDavis et all I1985T L Mer ger is thus an inevitable pro- 
cess and plays an essential role in galaxy formation 
and evolution. The potential impacts of mergers on 
galaxy evolution have been intensely discussed. In ad- 
ditional to being responsible for mass growth in the as- 
sembly history of galaxies, they a re, for example , ex- 
pected to t rigger quasar activities (ICarlberdll990t ) and 
starbursts (|Barnes fc HernquTstl 119911 ) . and thought to 
shape many important observational properties of galax- 
ies, such as star fo rmation rate, color and morphol- 
ogy transformation (|Toomre fc Toomrel fl972) . Galaxy 
mergers continue to be an important component for 
galaxy evolution from more recent works o n star for- 
mation history and remnant properties (e.g. ICox et ahl 
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| 2006f). starbursts. q uasars, and black hole growth (e.g. 
Hopk ins et al. | l2006l). bulge growth and morphology (e.g. 
i Hopkins et al.| l2010al) . baryonic content of galaxy (e.g. 
IStewart et al.ll2009af ) , and formation of dif fuse intra- halo 
light components (e.g. iPurcell et al.ll2007l ). 

Many studies have explored the galaxy merger rate 
as a function of cosmic time to understan d the assem- 
bly h i story of galaxie s obser yationally (iPatton et all 
12001 iConselice at al.| [2001 iLin et al.l 12004 1 2008] 

iKampczvk et al.l 12 007: 
iHsieh et al.l 120081: iBundv et atl 



Lotzet all l2008al [ 2011: 



Kartaltc pe et al.l 12 007: 
2009i jde Ravel et al 



2009, 120111: iChou et al.l 2011 
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Xu et al.l 120111) and theoretically (e.g. , iBerrier et al 
20061 lMateusTl2008T: IStewart et alj|2009bl: IHopkins et al 



l2010bf ). Since galaxy merger assembles two or more 
galaxies into one, it is expected to be more common 
in dense regions, i.e., having strong dependence on 
environments. To what degree this can be so i s an issue 
of gre at concern. The observational work of ILin et al.l 
(2010) was the first study to address this issue. They 
probed the environments of wet, dry, and mixed galaxy 
mergers at 0.75 < z < 1.2 using close pairs in the DEEP2 
Galaxy Redshift Survey, and found that the environment 
dependence of the galaxy merger rate for wet mergers 
is marginal while the dry and mixed merger rates 
increase rapidly with local density. A similar conclusion 
that merger rate has stron g environment depend ence 
latter was also obtained bv Ide Ravel et all (j 2 1 lh and 
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iKampczvk et al.l (|2011[) . However, galaxy merger rate 
in observation is actually an indirect physical quantity 
and can only be obtained by combining the fraction of 
galaxies that reveal signatures of merger activity from 
observations, and the merger time scale and the merger 
probability estimated from theoretical models. The 
environmental dependence of the merger probability and 
the merger timescale in theoretical models is less ad- 
dressed especially. Hence, to interpret the observational 
results more precisely, theoretical determination on how 
the merger time varies with environments and on how 
much pair contamination is due to projection effect in 
terms of environment is essential. There have been only 
several theoretical studies to address environmental 
dependence associated with mergers (iFakhouri fc Mai 
20091: [Hester fc Tastisiomill201d : iTonnesen fc Cenll2011h . 
Fakhouri fc Mai (|2009f ) utilized the Millennium Simula- 
tion to derive the merger rate of FOF DM halos as a 
function of the local mass density within a sphere of 7 
h^ 1 Mpc, and demonstrated that mergers occur about 
2.5 times more frequently in the densest regions than in 
voids at both z — and higher redshifts. However, their 
study focused on mergers of distinct FOF halos, while 
galaxy mergers do n o t nece ssarily follow their scenario. 
iHester fc Tastisiomil (|2010f ) alternatively adopted the 
merger tree of subhalos from the Millennium database 
to explore the same issues. The local halo environment 
in their study is defined as the count of all distinct halos 
and sub-halos, with V ma x < 120 kms^ 1 within a sphere 
of comoving radius 2 h^ 1 Mpc centered on each halo. 
They found that the merger rate is roughly independent 
of environment. However, the evolutionary history of 
subhalos may not exactly be that of the galaxies, due 
to different treatments on on orphan galaxies when 
subhalos are disrupted which may lead to different 
conclusions. The requirements for a relation between 
the histories of a given dark matter halo and of a galaxy 
inside that halo make the conversion from dark matter 
(halo/subh alo) merger rates to galaxy merger rates 
non-trivial (jHopkins et al.ll2010bft . 

In this study, we make use the publicly available galaxy 
catalogs based on four different semi-analytical mod- 
els (SAMs) implemented in the Millennium database. 
In addition to exploring the environmental dependence 
of galaxy merger rate, we also investigate the merger 
probability of close pairs (C mg ), which provides a cor- 
rection due to the projected selection criterion, and 
the merger timescale (T mg ), which provides the mean 
merger time for a close pair to merge, as a function 
of environments. These two quantities are important 
for observations to convert the pair fraction into galaxy 
merger rate. Besides, it is illuminating to understand 
how different SAMs affe ct galaxy merger rate a nd these 
two quantities. In fact, iHopkins et al.1 (|2010bft already 
pointed out that different theoretical methodologies lead 
to order-of-magnitude variations in predicted galaxy- 
galaxy merger rates. Whether such variations among 
different SAMs propagate into the environment depen- 
dence requires careful scrutiny. 

The paper is structured as follows. We briefly describe 
the Millennium simulation and the four galaxy catalogs 
from the semi-analytical models in use in Section 2. Sec- 
tion 3 is divided into five sub-sections. Environmental 
estimators are introduced in Section 3.1. In Section 3.2, 



galaxy merger rate is defined, where the redshift evolu- 
tion of the pair fraction N c and galaxy merger rate be- 
tween the models and observations are compared. The 
environmental dependence of N c and galaxy merger rate 
is presented in Section 3.3. Our results on C mgi and T mg 
as a function of environment are shown in Section 3.4. In 
Section 3.5, the impact of mass ratio of merging galaxies 
on the merger rate in different environments is discussed. 
Finally, discussion and conclusions are given in Section 
4 and 5, respectively. 

2. MODELS 
2.1. Simulations 

We make use of the Millennium Simulation, provided 
by the Virgo Consortium, in this study. The Millen- 
nium Simulation follows the hierarchical growth of dark 
matter (DM) structures with N = 2160 3 particles of 
mass 8.6 x 10 8 hr 1 M@ in a co-moving periodic box of 
size 500 hr 1 Mpc on a side fro m redshift z — 127 to 
the present (|Springel et al.ll2~005ft where the Hubble con- 
stant is parameterized as Hq — 100 h km s _1 Mpc -1 . 
The concordance ACDM cosmology is assumed in the 
Millennium simulation with parameters VL m = 0.25, 
Q b = 0.045, h = 0.73, Q A = 0.75, n = 1, and <j 8 = 0.9, 
consistent with a combined analysis of the 2dF Galaxy 
Redshift Survey and the first-year Wilkinson Microwave 
Anisotropy Pro be (WMAP) data. The GADGET-2 code 
(|Springe]| [2005ft . a variant of the TreePM method, is 
adopted for evaluating gravitational forces, allowing for a 
very large dynamic range and high computational speed 
even in situations where the clustering becomes strong. 
With a soften length of 5 h^ 1 kpc on a comoving scale, 
taken as the spatial resolution limit of the computation, 
the effective dynamic range thus reaches 10 5 in spatial 
scale. In addition, haloes hosting all luminous galax- 
ies brighter than 0.1 L„ (where L* is the characteristic 
luminosity of galaxies) can be resolved with a mass res- 
olution of ~ 100 particles. There are 64 epochs of out- 
put data spaced approximately logarithmically in time 
at early epoches and approximately linearly in time at 
late epoches (with A f ~ 300 Myr). 

A friends-of-friends (FOF) group finder ([Davis et alJ 
11985ft with a linking length of b — 0.2 in units of the 
mean particle separation is used to identify halos in the 
simulation. Each FOF identified halo is further broken 
into constituent subhalos (each with at least 20 parti- 
cles or 2.35 x 10 10 hr 1 M &) by the SUBFIND algo- 
rithm ([Springel et a l. 2001) which identifies gravitation- 
ally bound substructures within the host FOF halo; SUB- 
FIND typically finds a large background host subhalo 
and a number of smaller satellite subhalos. With all ha- 
los and subhalos determined, the hierarchical merging 
trees containing the details of how structures build up 
over cosmic time can then be constructed. These trees 
are the key information needed to compute the physical 
properties of the associated galaxy population for SAMs. 
It is known that the merger tree construction basically 
is to connect subhalos across all snapshot outputs, and 
established through the determination of a unique de- 
scendant from any given halo. In the Millennium sim- 
ulation, each potential descendant is scored based on a 
weighted count giving higher weight to particles that are 
more tightly bound in the halo under consideration, and 
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the unique descendant is the one with the highest score. 
In this way, the merger trees are defined and constitute 
the basic input needed by SAMs. 

The other simulation, ACDMiqqi,, adopted in this 
work is evolved in the fiat ACDM model: flo = 0.3, J7a 
= 0.7, D, b = 0.05, h = 0.7, and a s = 0.94, using 
512 3 DM particles in a co-moving periodic box of size 
100 h~ x M-pc on a side from redshift z = 100 to the 
present (|Jian e t al. 200§). This simulation also employs 
the GADGET-2 code, and the softening length e is set 
to 10 h^ 1 kpc (comoving) before redshift z — 2.3 and 
modified to 3 hr 1 kpc (physical) thereafter. Addition- 
ally, the outputs are at 67 epochs spaced equally in 100 
h^ 1 Mpc comoving length, and therefore, the time dif- 
ference for two adjacent snapshots increases slightly. For 
example, Af~ 0.24 Gyr at z ~ 1 and A t ~ 0.4 Gyr at 
z ~ 0. 

In the simulation ACDAfioob, a sub halo finder called 
HiFOF (Hierarchical Friends-of-Friends: IJian et al.ll2008| ) 
is hired to identify bound substructures. HiFOF first 
finds halos with a linking length of b = 0.2 in units of 
the mean particle separation (the lowest level), and starts 
to shorten the linking length linearly like a zoom-in pro- 
cess and perform virial condition check on those newly 
found substructures satisfying a criterion that the mini- 
mum particle number for subhalos has to be greater than 
30, corresponding to 1.9 x 10 10 hr 1 A/©. The process 
will continue for both those bound or unbound structures 
until it reaches a certain level or find no further bound 
structures. The merger tree in this simulation is defined 
differently from the Millennium one in how the unique 
descendant is determined. The criterion is that poten- 
tial descendants are the subhalos at least with 30 percent 
members overlapped with their progenitor members, and 
the descendant with most particles is the unique descen- 
dant. We also test the same tree definition as the one in 
Millennium and the result turns out to be similar. For 
example, the merger rate derived from these two trees dif- 
fers less than 4% for any redshift. Moreover, the galaxy 
model applie d to the simulati on ACDMiqqi is a simpli- 
fied one (see Uian et all 120081 for detail) that subhalos 
residing in the host halos more massive than a certain 
mass ~ 4 x 10 11 h^ 1 M© at z — are galaxy samples, 
and others are dark halos. This threshold is not based 
on galaxy luminosity but on the mass of the host halo 
in which galaxies reside. Due to the simple treatment 
for the galaxy properties, no orphan galaxies are consid- 
ered, and the result from this simulation included in this 
analysis is to provide an extreme reference case relative 
to SAMs. The result allows us to understand whether 
the merger rate derived from subhalo mergers (without 
orphan galaxies) may be significantly different from the 
merger rate derived from galaxy mergers. 

2.2. Semi- analytic models 

Semi-analytic modeling on galaxy formation and 
evolut ion was initially proposed by iWhite fc Frenkl 
(200 ll), and since then variations of the origi- 
nal SAM ha ve been proposed a nd verified to be 
a successful (|Croton et al.l I2006L iBower et al.l 120061: 
de Lucia fc Blaizotl 120071 : iBertone. De Lucia. &; Thomas! 
2007t iFont et al.l 120081 ). In this study, we focus on 
the galaxy catalogs produced by four different SAMs 
implemented on the Millennium simulation including 



Bower06 (IBower et al.ll200ll. Font08jFont et al.| [200l, 
peLucia06~ tide Lucia fc Blaizotl 120071) . and Bertone07 
(|Bertone. De Lucia, fc Thomas! I2007D . We investigate 
how much the results converge for the galaxy merger rate 
as a function of environment among different models. In 
addition, these galaxy catalogs are utilized to obtain in- 
formation such as the merger probability of close pairs 
and the merger time which are essential for observations 
to determine the galaxy merger rate. These catalogs are 
publicly available on the Millennium download sitqj. In 
the following, brief descriptions on different SAMs are 
given. For detail, the readers are refereed to the these 
papers. 

Bower06 is based on the Durham semi-analytic model 
of galaxy formation, GALFORM, and it implements a 
new treatment for the active galactic nuclei (AGN) en- 
ergy injection, determined by a self-regulating feedback 
loop that is assumed to quench cooling flows in massive 
haloes. They found that the feedback mechanism nat- 
urally produces a break at the bright end of the local 
luminosity function. Besides, satellite galaxies are as- 
sumed to lose all of their hot gases to the new parent 
halo upon being sufficiently close to the parent halos. 

DeLu cia06 is a slig h tly m odified version of the one 
used inlSpringel et al.l (120051) . iCroton et all (|2006l) . and 
Ide Lucia fc Blaizotl (120071) It in cludes the AGN feedback 
model from ICroton et al.l d2006) to suppress the cooling 
flows, and the model seems to be extremely efficient in 
switching off cooling in relatively massive haloes even 
at high redshifts. Additionally, they also employed the 
supernova feedback model from ICroton et al.l (|2006f ) to 
drive strong winds which blow out all gases away from 
galaxies on short time-scales. 

Bertone07 presents a new feedback scheme, which re- 
places empirical prescriptions on SN feedback of the Mu- 
nich semi-analytic model with a dynamical treatment 
of the galactic wind evolution, where ejection and re- 
cycling of gases and metals are treated self-consistently. 
It was shown that the observed mass-stellar metallicity 
and luminosity-gas metallicity relationships are able to 
be reproduced by their model. However, two drawbacks 
exist in their model. The number of bright galaxies in the 
LFs tends to be overestimated, and the color distribution 
of galaxies does not display the sharp color bimodality 
observed for galaxies in the local universe. 

Finally, in Font08 which is one of the Durham SAMs, 
the major change is that they merge a model of strip- 
ping processes based on detailed hydrodynamic simula- 
tions into GALFORM semi-analytical model for galaxy 
formation, in an attempt to fix the problem that satellite 
galaxies in groups and clusters are redder than observed. 
They found that the effect of the ram-pressure stripping 
on the colors of satellite galaxies is significant. With 
assumptions of the gradual stripping as opposed to the 
sudden stripping in Bower06, a considerable fraction of 
hot gases in satellite galaxies is preserve for several Gyr, 
and satellite galaxies remain blue for a relatively long 
p eriod of time as a resu lt. 

iKitzbichler fc White! (|2008t l noticed that different 
treatments in SAM influence the merger rate indirectly 
when identifying the merging systems as the observed 
merging galaxies. The galaxy formation modeling does 

5 http:/ /galaxy-catalogue. dur.ac.uk:8080/Millennium/ 
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a) at z- 1.0 b)atz-0.2 
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Fig. 1. — The distributions of the local density (l+5n) are plotted 
in different (the host halo mass) bins represented by color 
lines at two redshifts (a) z = 1.0 and (b) z = 0.2. It is observed 
that for a given overdensity range, galaxies are from various host 
environments, such as filed, group, or cluster environment. At high 
density region, the main contribution is from galaxies in group or 
cluster environment. On the contrary, filed galaxies contribute 
significant amount at low density area. It is also seen that the 
relative contributions from galaxies in different environment evolve 
with redshift. Two environmental estimators are found to be well 
correlated. 

not alter the dynamics of the underlying distribution of 
DM halos and subhalos, but when a galactic subhalo is 
tidally disrupted near the center of a more massive halo 
and eligible to merge into that halo, the model estimates 
the merger time based on a dynamical friction time ar- 
gument rather than assuming instant merging. In other 
words, once the subhalo disrupts, the galaxy evolution 
model waits for one dynamical friction time before the 
associated galaxy of this halo completes merging into the 
central galaxy of the main halo. After subhalo disrup- 
tion, the associated satellite galaxy is assumed to be at- 
tached to some particles, which are the most strongly 
bound particles of its last identified subhalo. The con- 
dition to obtain real i stic c lose pairs is demonstrated in 
IKitzbichler fc White! {2008) via comparisons on the two- 
point correlations between simulated galaxies and real 
galaxies. However, we notice that different treatments 
of orphan galaxies affect not only the time scales when 
a close pair completes merging but also the fraction of 
close pairs that will actually merge (see Section I3T4")) . 

3. GALAXY MERGER RATE AS A FUNCTION OF 
ENVIRONMENT 

3.1. Environmental estimators: Mh and (1 + 5 n ) 

To study the environmental effect, we adopt two dif- 
ferent environmental estimators. One is the host halo 
mass Mh for the merging galaxies, which is a clear phys- 
ical variable to represent the environmental effect and 
can be obtained easily from the simulation. However, 
the host halo mass is not a direct observable, and hence 
an environmental observable related to the local over- 
density density (1 + 5 n ) is also studied in conjunction 
with Mh- The local galaxy density (1 + S n ) is esti- 
mated using the projected n" l -nearest-neighbor surface 
density, £„ = n/irDp n , where D Pi . n is the projected 

distance to the n^-nearest neighbor that is within the 



line-of-sight velocity interval of 1000 km s _1 . In addi- 
tion, each density value is divided by the median S n 
of galaxies at that redshift over the whole sample, i.e., 
1 + S n = £„ / ' median{T, n ) to reduce the inf luence of vari- 
ation s in redshift sampling in the survey (jCooper et al.1 
l2005h . Note that, (l + 5 n ) is a relative but not a absolute 
quantity, and direct comparison for the same (1 + S n ) at 
two different redshifts is, thus, not appropriate. In or- 
der to make dire ct comparison with DEEP2 results from 
iLin et ail (|2010| ). n = 6 is chosen to account for using 
the 3 rd -nearest neighbor in the DEEP2 sample with the 
average redshift completeness ~ 50%. In the rest of this 
paper, we refer 1+S n to the overdensity measured up 
to the 6 t?l -nearest-neighbor for galaxies. Furthermore, 
the applied galaxy selection cuts are different in the sim- 
ulation ACDM wob and in the SAMs. For the SAMs, 
galaxy samples are selected with a flux limit cut R < 24.1 
which is the same cut as in DEEP galaxies. For simula- 
tion ACDMioob, the median distance to the n'^-nearest- 
neighbor above certain mass cut M cut for subhalos is 
computed so as to match the median distance to the 3 rd - 
nearest neighbor in the DEEP2 sample. The resulting 
(1 + S n ) distribution of halos for ACDMioob displays a 
very similar profile with the (1 + S3) di stribution of ob - 
served DEEP2 galaxies (see Figure 4 of ILin et al.ll2010lh 
The corresponding M cut for n = 6 are (a) 2.23, (b) 2.48, 

(c) 5.75, and (d) 8.35 x 10 10 fr" 1 M for four different 
redshifts, i.e. z = (a) ~ 1.1, (b) ~ 0.9, (c) ~ 0.7, and 

(d) ~ 0.5. 

With these two estimators Mh and (1 + S n ), it is 
interesting to understand how they are correlated and 
how their correlation evolves. Therefore, the local den- 
sity (1 + 5 n ) distributions for galaxies resided in hosts 
of different masses denoted in color-coded lines from 
Font08 are plotted in Figure [T] as an example. We 
split a total sample of N galaxies at a certain redshift 
into four bins based on the host halo mass of galax- 
ies. These M h bins are 11 log w (M h /( hr 1 M Q )) < 

12 (red), 12 «C log w (M h /( h' 1 AT©)) < 13 (green), 

13 log w (M h /{ h" 1 M )) ^ 14 (blue), and 14 sC 
logio(Mh/ ( hr 1 M Q )) ^ 15 (black). In each of these 
four Mh bins, we further divide galaxies into different 
(1 + 8 n ) bins. That is, the total sum of galaxy counts 
from all (1 + 8 n ) bin and from four Mh bins should equal 
to N. We subsequently use N as a normalization factor so 
that the sum of normalized count from all (1 + 8 n ) bins 
for the light blue line is equal to 1. It is noted that for a 
given range of local density (1 + S n ), its total strength is 
a collective result combing different contributions from 
various hosting environments. In other words, a galaxies 
with the highest (1 + S n ) is not necessary from the most 
massive hosting environment. Nevertheless, it is seen in 
Figure [1] that for those galaxies in more massive host 
halo, their local density (1 + 8 n ) tends to be distributed 
in a denser regions. In addition, it is observed that when 
comparing the local density distributions at two differ- 
ent redshifts that the relative contribution from different 
Mh bins evolves with redshift. The contribution from 
the two most massive Mh bins increases with z and that 
from the two least massive bins decreases with z, im- 
plying that in time field galaxies gradually merge with 
galaxy clusters. Therefore, we conclude that two estima- 
tors are well-correlated and the local density indicator 
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(1 + 5 n ) has a real physical meaning for the environment, 
and the dominant contribution for large (1 + S n ) is from 
groups or small clusters. 

3.2. Pair fraction N c and galaxy merger rate 

There are two observational approaches used to probe 
the evolution of merge rs. One is relat e d to the 
close-pair c ount (e.g. iPatton et al.l 120021 iLin et all 
12004 120081) . and the other is to count galaxy mergers 
throu gh morphologica l signatures of galaxy in teraction 
(e.g. IConselice at all 120031 ; ILotz et all I2008al) . In this 
study, we follow the close-pair count approach to study 
the galaxy merger and explore its environmental depen- 
dence. In this approach, the direct observable is the av- 
erage number of companions per galaxy, defined as 



N, 



_ 2N pai 



(1) 



where N pa i r is the number of galaxy pairs and N g is 
the number of galaxies in the sample. Through study- 
ing N c from different models, a direct comparison be- 
tween theories and observations can be made. Normally, 
the pair fraction is not equal to the average number of 
companions per galaxy, N c , but in this study, we liber- 
ally refer the term "pair fraction" as "the average num- 
ber of companions per galaxy" . Follow ing the observa- 
tional criterion in IPatton et al.1 ()2002[ ). projected close 
pairs are so defined that the projected separation satis- 
fies 10 hr 1 kpc < Ar < r max , where r max — 50 h~ x kpc 
is assumed in this study, and the rest-frame relative ve- 
locity Av is less than 500 km s _1 . Additionally, galaxy 
samples are selected in the evolution-corrected B band 
magnitude range,— 21 < Mg < —19, where Mg is de- 
fined as Mr(z = ) + Qz, with Q = 1.3 adopted from 
iFaber et all |2007T ). 

Our results on the evolution of the pair fraction N c 
are shown on the left panel in Figure [2] This fig- 
ure contains observa tion results (red erro rbars) taken 
from lLin etafl (12004) and ILin et"atl ([20081 ) that include 
DEEP2 (jDavis et al.ll2003L 120071 ) and some low redshift 
data, and simulation results from four semi-analytical 
models, Bower06, Font 08, Delucia06, a nd Bertone07, 
alon g with results from Uian et al.l (I2008D (AC DM WQb ) 



and iBerrier et ahl (|2006l ). iBerrier et I all (j200ll adopted 
two simple models for associating subhalos with galaxies: 
Vi n , the maximum circular velocity that the subhalo had 
when it was first accreted into the host halo, and V now , 
the maximum circular velocity tha t the subhalo has at 
the current epoch. Two models of IBerrier et al.l (2006) 
give a reasonable range compared with observations, as 
shown in Figure [2J It appears that the theoretical N c (z) 
vary with models by a factor as large as ~ 3 and over- 
all slightly deviate from the observations. However, the 
flat evolutionary trends in the models are consistent with 
the observations. The flat trend of N c (z) was first the- 
oretically obtained by IBerrier et all (|2006l) and our re- 
sults from the models also reveal similar trends, suggest- 
ing that both observations and simulations support weak 
redshift dependence in N c . Additionally, the projected 
close-pair counts are related to the projected two-point 
correlations, and hence stronger clustering implies more 
pair counts. To further examine the deviations found in 
N c , the clustering of the galaxies at small scale satisfy- 



ing the magnitude selection criterion from the models are 
explored. The projected two-point correlation functions 
w p (r) of different models are computed for comparisons 
and shown on the right panel in Figure [2j It is seen that 
the clustering strength at small length scale varies with 
different models, and Bower06 and Font08 show stronger 
clustering than DeLucia06 and Bertone07 at length scale 
r p = 50 h^ 1 kpc. The higher strength found at small 
length scale in Bower06 and Font08 indicates a higher 
value of N c in these two models. 

From observational point of view, galaxy merger rate 
is not an observable but an indirect quantity converted 
from the pair fraction N c . On the contrary, the theoret- 
ical merger rate is a direct measurable in simulation. By 
taking this advantage, we connect the theoretical merger 
rate to the merger rate in observational form to derive 
T rng and to know the real meaning of T mg defined in ob- 
servation. We start from giving a clear definition on how 
the merger rate is measured in this work. There are two 
definitions considered for the galaxy merger rate. The 
first one is the normalized merger rate, and it is defined 
and estimated as follows. 

N„ 



Fn 



2 x 



i=i 1 
N„ 



(2) 



where N mg is the number of mergers, Ti is the time to 
form the i th merger when the merger was still an iden- 
tified close pair. In Equation [2] the factor of 2 is intro- 
duced to turn the number of merger pairs into the num- 
ber of merging galaxies. From Equation [21 F mg shows 
dependence on the Ti distribution of the merger pairs. 
To quantify the Tj distribution, we follow the identified 
close pairs at four starting redshifts until they merge and 
record its distribution. The results are presented in Fig- 
ure[3l The distributions among different SAMs are nearly 
flat and agree with each other. The flat distribution re- 
sults in a longer medium or mean merging time about 3 
to 4 Gyrs. On the contrary, the distribution from the 
simulation ACDMiqqi, shows a peak at short merging 
time. The main reason for the deviation is correlated 
with the lack of orphan galaxies in ACDM100& and the 
merging time turns out to be much shorter than in the 
SAMs. 

In contrast, the definition from lLin et all (|2010l) for the 
normalized merger rate is 



f mg = {l + G) CmsNc 



T, 



(3) 



where G is the correction factor that accounts for the 
selection effect of companions due to the restricted lu- 
minosity range, and C mg represents the probability of 
galaxies in projected close pairs that will merge before 
the present. In mathematical expression, 



Gmg J Npav 



(4) 



Because not all projected pairs are merging systems, C mff 
is introduced to account for the contamination from in- 
terlopers due to the difficulty in disentangling the Hub- 
ble expansion and the galaxy peculiar velocity. To obtain 
fmg, Cmg is hence as important as T mg , and these two 
essential quantities are not direct observables but can be 
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evaluated through simulations. Apart from the factor of 
(1 + G), F mg and j mg should be the same. The normal- 
ized merger rate F mg can then be expressed in terms of 
N c , C m g, and T mg . Substituting C mg with Equation |4] 
and N c with Equation [TJ we then obtain that 



F, 



mg 



1 m.n 



2N, 



mg 



N T 



(5) 



Rewriting Equation [5J T mg can then be expressed as a 



function of F rng , N g , and N mg such that 



T 



2N„ 



i—l 



(6) 



The determination of T mg is based on this formula in our 
work. 

There is an innate difference between the true merger 
timescale, and the timescale over which the merger would 
be obser vable, i.e. the "ob servability timescale" , see for 
example lLotz et all (|2011| ). T mg is also an average ob- 
servability tnnescale, butdefined differently from Equa- 
tion (8) in lLotz et al.l (1201 lh . Our T m „ definitio n is close 
to Equation (8) in iKitzbichler fc Whhel (|2008[ ). and the 
main difference is that they absorb C mg into the merger 



time but we do not, making their merger timescale longer 
than ours. 

The other definition is the volumetric merger rate 
which has a form that 



I? 



2 — 1 



mg 



L 3 



ljF mg x rig(z), 



(7) 



where L is the length of the simulation box on a side 
in unit of h^ 1 Mpc and n g (z) is the comoving number 



density of galaxie s, i.e. 



N g /L 3 in this study. 



Comparatively, in iLin et al. (2004|) it gives that 



T mn = (0.5 + G) x n v C ^ Nc 



mg 



T, 



(8) 



'».</ 



where G is the correction factor previously described and 
the factor of 0.5 in (0.5+G) converts the number of merg- 
ing galaxies into the number of merger events. Except 
for the factor of (0.5 + G), R mg is the same as T mg . 

In Figure 31 the normalized merger rate F mg (top) 
and the volumetric merger rate R m g (bottom left) for 
the models and observations are plotted as a function 
of z, plus the evolution of the number density of the 
selected galaxies n g (bottom right). Deviations in the 
merger rates among the models are seen, but the flat 
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Fig. 3. — Distribution of the merging times in Gyr at four different redshifts. Four SAMs are in good agreement. The distribution of the 
merging times gives a approximately constant profile from the SAMs. while there is a peak in the profile from ACDMiqqi,. Galaxy close 
pairs in AC DAIiqqi, have a large fraction of pairs that merge quickly leading to a peak at short merging time. The reason for the short 
merging time is mainly due to the lack of treatment on the orphan galaxies. 



evolutionary trends in the merger rate for the models 
are consistent in the figures of F mg and Rmg- With the 
parametrization F mg (z) oc (1 + z) m in rcdshift range 
0.2 < z < 1.0, we find that m ~ 0.71 ± 0.24 for Font08, 
- 1.88 ±0.19 for Bower06, - 0.359 ±0.31 for DeLucia06, 
and ~ 0.0±0.28 for Betone07. When R mg is also assumed 
to be oc (1 + z) m , we find m ~ 0.46 ± 0.38 for Font08, 
~ 0.93 ± 0.34 for Bower06, - 0.21 ± 0.33 for DeLucia06, 
and ~ 0.0 ±0.34 for Betone07. Besides, the large discrep- 
ancies between the models and observations are found 
mainly due to value s of C mn and T m „ adopted by observa- 
tions. For example. iLin et aLl (|2004l ) adopts a normaliza- 
tion constant C m g/T mg = 0.5/0.5 = 1 while it is ~ 0.15 
from Font08 and Bower06, ~ 0.2 from DeLucia06, and 
^0.17 from Bertone07. This also reflects a fact that the 
merger rate in the semi-analytical models implemented 
in Millennium simulation is comparatively smaller than 
those in other simulations. Moreover, the mass ratio of 
the merging galaxies has impact on the galaxy merger 
rate, with much higher merger rates when smaller galaxy 
mass ratios are considered, making comparison from ob- 
servations to theory very uncertain. On the other hand, 
the result from the simulation AC D Mioob appears to be 
broadly consistent with observational data points with 
slopes m ~ 2.18 ± 0.11 in F mg and ~ 1.72 ± 0.2 in 
R m g- The consistency arises primarily from the nor- 
malization constant C mg /T mg found in ACDAfioob to be 
close to that adopted in observations. The simulation 
ACDMiQQb uses a simplified galaxy model without con- 
sidering the orphan galaxies that have already lost their 



surrounding subhalo, and thus it lacks contribution from 
these orphan galaxies when estimating the merger rate. 
It should be noted that the treatment for these orphan 
galaxies in semi-analytical models is such that once the 
subhalo disrupts, the orphan galaxy waits for a dynami- 
cal friction time before merging into the central ga laxy of 
the main halo (e.g. seelKitzbichf er fc Whitdl2008l ) . That 
is, the merger timescale in the SAMs is much longer than 
that in ACDMi 00 b- Consequently, the merging history 
turns out to be very different, yielding in a higher merger 
rate in AC-D-Minnft. Howe ver, it was demonstrated by 
iKitzbichler fc White! ((2008) that with the orphan galaxy 
treatment, clustering at small scale at z = produces 
a similar correlation strength as in SDSS whereas there 
will be a deficit in the correlation strength without the 
treatment. The consistency of ACDMioqi, may therefore 
seem to be a coincidence and subhalo-subhalo mergers 
may present only a partial picture of real galaxy merg- 
ers. Our inclusion of the results from AC-DM100& in this 
study is simply to provide a reference. 

On the bottom-right panel in Figure HJ the evolution 
of n g with galaxies selected in two magnitude range is 
plotted for the models. The densities among the models 
can differ by an order of magnitude at high redshift and 
converge at low redshift. This plot is to demonstrate that 
under the same selection criterions, density variations 
exist among different models. 



3.3. N c and F mg as a function of environment 
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basically is F„ 



In the previous section, redshift evolution of the pair 
fraction and the merger rate is discussed without con- 
sidering environmental effect. To assess the effect of the 
environmental dependence of the pair fraction, N c is plot- 
ted as a function of local density (1 + S n ) (left) and Mh 
(right) in four redshift ranges (a) 1.0 < z < 1.4, 
(b) 0.8 < z < 1.0, (c) 0.5 < z < 0.8, and (d) 
0.2 < z < 0.5 in Figure [5] When we bin galaxy close 
pairs into different environments, we exclude pairs whose 
pair galaxies are in different (1 + S n ) or Mh bins, and we 
find that the fraction of the eliminated pairs in Font08 
for Mh bins, for example, is ~ 15% at z ~ 1 and drops 
to ~ 6% at z = 0. The fraction is small and unlikely to 
affect the result. On the left panel, it is seen that the 
pair fraction iV c increases with the local density (1 + S n ) 
and different models are in g ood agreement with DEEP2 
observation (|Lin et al.l|2008f ), except for the densest bin. 
The increasing tendency is understandable simply be- 
cause galaxies in overdense environment are more clus- 
tered, and thus easily appear in pairs. By contrast, on 
the right panel N c appears to grow at low Mh bins, but 
gradually gets flat and displays large deviations among 
the models at high Mh bins. The explanation for this is 
that any given bin of (l + 8 n )) is contributed from various 



Mh bins, as discussed in Section HO (see Figure [T]), and 
due to least contribution from the most massive cluster 
scale in overdense regions (the bins with high (1 + <5„)), 
the total effect leads to a result that N c continues to grow 
as local density (l + <5„) increases. Additionally, the evo- 
lution of N c appears to evolve weakly with redshift z in 
any case. However, different SAMs differ considerably 
by a factor of ~ 1.5 in N c in high density region and by 
a factor of ~ 4 in the highest Mh bin. 

The normalized galaxy merger rate F mg is also evalu- 
ated in terms of (1 + 5 n ) (left) and Mh (right) in Fig- 
ure [6] On the l eft panel, the observational results from 
ILin et al.1 poll are included for comparisons. It is evi- 
dent that F mg has strong dependence on their surround- 
ing environment, yielding that mergers occur more fre- 
quently in more dense region than in under-dense re- 
gion. In addition, F mg appears to be nearly redshift 
independent but model dependent. On the left panel, 
if the galaxy merger rate is parameterized as F mg oc 
(1 + S n ) a , it is seen that F mg is flatter in Bower06 (cyan) 
and Font08 (green) with a ~ 0.3, and steeper in DeLu- 
cia06 (pink) and Bertone07 (blue) with a ~ 0.6 to 0.7. 
That is, for all redshifts, galaxies merge more rapidly 
in high density regions than in under-dense regions by 
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a factor of ~ 4 in Bower06 and Font08 and by a fac- 
tor of ~ 20 in DeLucia06 and Bertone07. The posi- 
tive correlation between galaxy merger rate and the lo- 
cale densit y is in broad agreeme nt with the measure- 
ment from iFakhouri fc Mai (|2009( ) , in that for galaxy- 
mass halos, mergers occur ~ 2.5 times more frequently 
in the densest regions. The strong environment depen- 
dence is also c onsistent w i th th e results in recen t obser - 
vati onal works iLin et all (1201(1) . Ide Ravel etaLI (|2011l ). 
and lKampczvk et al.l (|2011| ) over the redshift range 0.2 
< z < 1.2. In contrast, F mg in ACDM 10 ob (black) shows 
nearly fiat profiles implying that it is roughly indepen- 
dent of environment. This finding of no environmenta l 
dependence is close to what iHester fc Tastisiomil ([2010D 
concluded using a merger tree of subhalos from the Mil- 
lennium simulation. The difference between the SAMs 
and ACDMioob will be further discussed in Section |U 

On the other hand, when F mg is expressed in terms of 
Mh, the profile shows a turn-over and has a peak in the 
Mh range between 10 12 and 10 13 hr x M Q , which cor- 
responds to group environments. The major disparity 
among the models is at the most massive Mh bin. Merg- 
ers occur approximately an order of magnitude more fre- 
quently in group or cluster environments than in field en- 
vironment, and the peak indicates that galaxies in group 
environments merge most efficiently. Because galaxies in 
a local density (l + 5 n ) bin are contributed partially from 
field, group, and cluster environments (see Figure Q}, the 
actual enhancement in the merger rate seen in high local 



density regions on the left panel of Figure [5] is mainly con- 
tributed from galaxies in group envir onment. Th i s find - 
ing is in good agreement with that of iTran et al.l (2008) 
that the group environment is critical for mergers to form 
massive galaxies. The fact that group galaxies have a 
high merger rate is expected, since in such environments 
galaxy close pairs have low velocity dispersion and are 
hence easier to merge. In addition, as we will discuss in 
Section 13.41 close pairs in group environments are also 
less contaminated from the project effect. 

3.4. Environment dependence of the merger timescale 
T mg and the merger probability of close pairs C mg 

Because C mg and T mg are two important quantities for 
observations to convert N c into galaxy merger rate and 
can be determined only theoretically, we explore them in 
detail in this section. We evaluate T mg from Equation [6] 
previously defined. In Figure O T mg is expressed in term 
of two environmental estimators, (1 4- 5 n ) (left) and Mh 
(right). On the both panels, the merger time T mg shows 
weak dependence on environments. In the densest re- 
gion, T mg is ~ 10% shorter than that in under dense 
region, and the deviation of T mg among different models 
is small. Moreover, the merger timescale T mg is much 
longer in the semi-analytical models than in AC DAI wot,. 
The treatment on orphan galaxies to add dynamical fric- 
tion time in the semi-analytical models is responsible for 
the longer timescale. We also observe that T mg depends 
slightly on redshift z and declines as (1 + z)^ 1 for SAMs. 
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Left: Frag increases with (1 + 5 n ), and in DeLucia06 (pink) and Bertone07 (blue) the increasing strength 



relatively flattened and is ~ 4 time larger. Right: F mg reveals a turnover profile, and the difference between that of field galaxies and that 
of group galaxies is over an order of magnitude. Galaxies in group environment appear to have highest merger rate. The high merger rate 
in high local density regions is attributed to the significant contributions from galaxies in group environment. 



However, we argue that the rcdshift dependence may not 
be a real effect. The shorter T mg obtained at low redshift 
is simply because many of projected close pairs identified 
at low z do have no time to merge before the present. 
That is, the shorter T mg represents only a small merger 
population. 

Regarding the evaluation on C mg , we follow the merger 
tree forward in time in the simulations to examine 
whether projected close pairs will merge before the 
present. If close pairs do not merger before the present 
time, they will not be counted as mergers. The merger 
probability of close pairs C mg is the fraction of the merger 
number to the close pair number and is to account for 
the effect of interlope contamination. 

The merger probability of close pairs C mg is plotted 
as a function of local density (1 + S n ) (left) and Mh 
(right) in Figure [5] Contrary to weak dependence on 
environments in T mg , C mg reveals strong dependence on 
environments. On the left panel, C mg monotonically de- 
clines with the local density (1 + S n ) while it has a peak 
at the mid M h bin (12 < log 10 (M h /M Q ) < 13), cor- 
responding to the group environment. The profiles of 
C mg in different SAMs also show significant discrepan- 
cies by a factor of 3 on both left and right panels. At 
low redshift, the C mg profiles become flatter; this is un- 
derstandable since a number of projected close pairs do 
not have enough time to merger before the present, and 
C m g does not vary with environment richness sensitively, 
thereby yielding a flatter profile and have lesser difference 



between underdense and overdense regions. In contrast, 
C mg at low z in ACDMioob still shows a steep profile 
and differs with those of SAMs by a factor of ~ 2 at low 
density regions. It is mainly due to the short merging 
timescale T mg in ACDMi 00 b (see Figure [3]), and result- 
ing in a narrow distribution of merging time, and hence 
there is little difference between high and low redshifts. 

To understand the origin leading to the environmen- 
tal dependence found in C mg , the ratio of N^£)/N va i r is 
investigate to give an idea of how many projected close 
pairs are really close in three dimensional space, where 
N3D is the number of galaxy pairs that are close in three 
dimensional space in the 2D close pairs sample. The ra- 
tio is computed as a function of local density (1 + 5 n ) and 
Mh, and the results are shown in Figure [51 The ratio of 
A3D /N pair displays curves similar to the profiles of C mg 
on both right and left panels. It is thus inferred that the 
projection effect is very likely the main cause responsible 
for the drop of C mg in overdense region. However, the 
slope in the profile of the ratio of N$d /N pa i r is not as 
steep as in that of C mg . It is possible that in the over- 
dense region projected close pairs are affected by some 
other physical processes, such as high velocity dispersion, 
to make merger harder to occur other than the projection 
effect. 

To investigate this issue, we also plot three-dimensional 
separation versus three-dimensional velocity difference 
for those 2D selected close pairs that really merge (top) 
and that do not merge before the present (bottom) 
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Fig. 7. — Merger time T mg as a function of local density (1 + 8 n ) (left) and (right). T mg is estimated using Equation [6] Left: T mg 



is nearly constant without environment (1 ■ 
or less shorter than in the low density region. 



At lower rcdshift, T„ 



L mg is apparently shorter. However, it doest not mean that projected 

will 

be obtained at low redshift. Right: T mg is also approximately independent of environment and the models. At high redshift, T mg in 
less massive bin (field environment) is twice longer than that in the most massive bin (cluster environment), but the errorbar in 



in Figure [TU] at z ~ 1.0, using Font08 as an exam- 
ple to know how galaxies populate in this two dimen- 
sional parameter space. In this example, the merger 
and non- merger probabilities of close pairs are ~ 31% 
and 69%, respectively. We find that for those merged 
pairs (~ 31%), about 78% of them satisfy the close-pair 
selection criterions that the physical separation is less 
than -y/3/2 50 h~ x kpc and the velocity difference is less 
than 500 km/s, and 22% of them are either with larger 
separation or large velocity difference but still merge in 
the end. On the other hand, for these non-merger pairs 
(~ 69%), about 36% of them satisfy the close-pair selec- 
tion criterions but do not merger in the end, and 64% of 
them that are not thrcc-dimensionally close and do not 
merge which can be attributed to the projection effect. 
The results suggest that the projection cause is the main 
cause for the close-pair contamination. However, sur- 
prisingly there are 25%(= 69% x 36%) 3D-close pairs do 
not merge, suggesting that it is indeed high relative ve- 
locities between the 3D-close pairs to prevent them from 
merging. 

3.5. the impact of mass ratio of merging galaxies on the 
merger rate in different environments 

There have been many discu ssions of merger rates 
and merger mass r atios (e.g. iFakhouri fc Mai 120081 : 
iStewart et aT]|2009bl: iHopkins et al.ll2010bD . To under- 
stand the impact of mass ratio of merging galaxies on the 



merger rate in different environments due to the luminos- 
ity cut sample selection in our study, we first estimate the 
fraction of pair satisfying the mass ratio criterion to total 
pairs as a function of environment in our luminosity cut, 
and we also examine its dependence on the mass ratio 
criterion. When we refer to merger mass ratios for pairs, 
we use the definition that \l = (Mg/Mg), the stellar mass 
of galaxy 2 over the stellar mass of galaxy 1, always taken 
such that < \i < 1. In the top panel of Figure [Til the 
fraction of merger mass ratios satisfying /i > 1/2, 1/3, 
1/4, and 1/6 to all merger pairs is plotted as a function 
of Mh ■ The fraction exhibits roughly no environment de- 
pendence for the four merger mass ratios, implying that 
there is no significant discrepancy for the distribution 
of /i across different environments. Moreover, the frac- 
tion also shows no redshift evolution, suggesting that the 
passive evolutionary luminosity cut should select pairs of 
different epochs with similar physical properties. In ad- 
dition, the galaxy merger rate F mg as a function of Mh 
for different merger mass ratios is plotted in the bottom 
panel of Figure [TTJ We find that F mg increases as fx de- 
creases, but the shape of F mg does not change. The sim- 
ilarity among different F mg profiles again indicates that 
the impact of mass ratio on the merger rate in different 
environments should be a minor effect. 
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Fig. 8. — Merger fraction C m g as a function of local density (1 + <5 n ) (left) and (right). C mg is not only affected by the density, but 
also by the models. Left: Cmg declines with the density. The variation between the lowest density and the highest density region is able 
to reach a factor of ~ 5 in Bower06 and Font08, but a factor of ~ 1.5 in DeLucia06 and Bertone07. At the same density bin, different 
model deviates as large as a factor of ~ 3.5. Right: The highest value of C m g is in the bin with 12.0 < logio(Mh/M<?)) < 13.0 (group 
environment). The deviation among different models is as large as a factor of ~ 10. 



By using galaxy catalogs based on four SAMs, we find 
that the normalized galaxy merger rate F mg reveal strong 
environmental dependence and such dependence is sim- 
ilar at all redshifts but strongly model dependent. We 
observe that when the galaxy merger rate is expressed in 
terms of (l + 5 n ) as F mg oc (l + 6 n ) a , a ~ 0.3 in Bower06 
and Font08, and a ~ 0.6 — 0.7 in DeLucia06 and 
Bertone07. The merger rate F mg between low density 
and high density regions can differ by a factor of ~ 4 in 
Bower06 and Font08, and a factor of ~ 20 in DeLucia06 
and Bertone07. In contrast, F rng in ACDMioob, rep- 
resenting a case of subhalo-subhalo mergers, shows less 
or no environmental dependence. Despite of large dis- 
crepancy in environmental dependence of F mg among the 
SAMs, the incre asing tend is broadly co nsistent with the 
result found by IFakhouri fc Mai (|2009f ) that for galaxy- 
mass halos, mergers occur ~ 2.5 times more frequently 
in the densest regi ons, whereas the resu l t in A CDMiqq}, 
is similar to what iHester k, Tastisiomil (|2010D obtained 
that the merger rate does not correlate with environment. 
However, it may not be fair to ma ke directly comparisons 
between our results and that of IFakhouri &: Ma (2009) 
simply because our sample contains different halo popu- 
lations whereas they did not. The fact that the positive 
correlation between F mg and the local density is observed 
in the SAMs and no correlation in ACZM/ioofc arises 
mainly from competition between increasing N c and de- 
creasing C mg with the local density with T mg remaining 
roughly constant. The pair fraction N c in the SAMs has 



a steeper increasing slope than that in ACDMioob, while 
the merger probability of close pairs C mg in the models 
has a flatter decreasing slope than that in ACDMioo&- 
(We note that the enhanced N c at high local density in 
the SAMs is attributed to the presence of orphan galax- 
ies.) In the end, the combined effect exhibiting in F mg 
turns out to be a flat profile in KC DMy^b and a strongly 
increa sing profile in the SAMs. As iKitzbichler fc White! 
(2008) pointed out, inclusion of orphan galaxies is cru- 
cial when calibrating the conversion from pair counts to 
merger rates. 

Deviation in the C mg and T mg determination among 
different SAMs is another interesting topic to be dis- 
cussed. We notice that C mg and T mg adopted by ob- 
servations vary with the referenced simulations, and the 
four SAMs we analyze give smaller C m , „ and longer 
Tm.„ than those taken by ob servations (e.g.. lPatton et afl 
I2002t iLin et al.l 12004 [20081 etc). It leads to significant 
deviation in the merger rate between the observation 
inference and the model predictions. In our analysis, 
C mg ranges from 0.27 to 0.55 depending on the models 
and T mg is about 2 Gyr nearly no model dependence 
at z ~ l. However, T m o is assumed to be 0.5 Gyr 
in the observation work in iPatton et al.l (|2002l ). In ad- 
dition to the larger projection separation 50 hr x kpc 
taken in our study to derive T mg as opposed to 20 or 
30 h~ l kpc taken by observations, the merger timescale 
is actually longer in the SAMs compared wi th simula- 
tion results from what lLotz et al.l (2008b, 2010) obtained, 
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Fig. 9. — Ratio of N30 /N pa i r as a function of local density (1 + S n ) (left) and Mh (right). From this plot we are able to know in these 
2D selected close pairs how many of them are three dimensionally close. It is found that the ratio has environmental dependence similar to 
Cmg, and it is thus inferred that the projection effect is responsible for the environmental dependence of Cmg- However, the slope of the 
ratio is flatter than that of C m g, and it is suspected that some other physical mechanisms in dense environment prevent close pairs from 
merging. 



~ 1 Gyr. Our merging time distribution is wide and 
nearly flat (see Figure [3]), and its median is ~ 3 Gyr 
at z ~ 1. This value may be considered to be a rep- 
resentative merger timescale. If th e merger timescale is 
estima ted following the definition in Kitzbichler & White 
(2008) that (T merge ) = T mg /C rng , it can resul t in an even 
longer timescale ~ 2 to 4 Gyr. Additionally, lLotz et all 
( 2008b) gave a morphological analysis from a large suite 
of N-body/hydrodynamical gas- rich disc galaxy merger 
simulations processed through realistic radiative transfer 
models. For the same projected separation and velocity 
selection as in this paper, the time scale they obtained 
rang es from ~ 0.3 to 1.4 Gyr and the median is ~ 1 
Gyr. lLotz et ail ()2010l ) further explored the effect of gas 
fraction, and mergers with different gas ratio result in 
time scale ~ 0.9 Gyr for 1:1 baryonic mergers, and ~ 
1.2 Gyr for 3:1 baryonic mass ratio mergers. Therefore, 
depending on what simulation is adopted and what the 
definition for the timescale is, the uncertainty of merger 
rate can be as large as ~ 4. 

In addition to the uncertainty of the merger time T mg 
as discussed above, the merger probability of close pairs 
C mg also shows large variation depending on how they 
are calculated. Especially, the environmental depen- 
dence of C m a is addre s sed to a much lesser extent in 
the past. In iLin et al.l (|2010D . the decrease of C mg in 
high and low local density regions is by a factor of ~ 
4, approximately independent of z. However, this re- 
sult corresponds to the case of subhalo-subhalo mergers. 



From our analysis, C mg is apparently model and redshift 
dependent. The decline is by a factor of ~ 4 in Bower06 
and Font08 and ~ 1.5 in DeLucia06 and Bertone07 at 
redshift z ~ 1 while the decline is by a factor of ~ 2 
in Bower06 and Font08 and nearly flat in DeLucia06 and 
Bertone07 at z ~ 0.4 . The uncertainty thus introduced 
by C mg could be as large as ~ 2 to 3. 

5. SUMMARY 

We make use of publicly available galaxy catalogs from 
four SAMs, including Bower06, DeLucia06, Bertone07, 
and Font08, implemented in the Millennium Simulation 
to explore the relation between the galaxy merger rate 
and its underlying environment, as well as to evaluate 
the model dependence of the merger rate. The approach 
taken by us closely follows the close pair observations 
by giving selections with an evolution-corrected B band 
magnitude range, —21 < Afjj < —19, with the projec- 
tion separation in physical length Ar between 10 and 
50 hr x kpc, and with the line-of-sight velocity difference 
less than 500 km/s. In this study, environment is quan- 
tified using the local density (1 + S n ) and the host halo 
mass Mh- Given these two estimators, we have been able 
to investigate the impact of environment on the galaxy 
merger rate. The results are summarized as follows. 

(1) The whole (1 + S n ) distribution is decomposed into 
four bins based on the host halo mass Mh that galaxies 
reside in to understand how these two environment es- 
timators are correlated and how the correlation evolves 
with redshift. It is found that galaxies in more mas- 
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Fig. 10. — An illustration on how galaxies populate in the two 
dimensional parameters space of 3D physical separation and 3D 
velocity difference for those 2D selected close pairs that merge (top) 
and that do not merger before the present (bottom) at z ~ 1 using 
Font08 as an example. On the top panel, the color bar indicates 
the merging time of the close pairs. Additionally, on top-right 
corner the red-color numbers 31% and 69% indicate the merger 
and non-merger probability of close pairs, respectively. The blue- 
color numbers inside and outside the box show the ratios that pairs 
satisfy and do not satisfy the projected close pair criterion. It is 
found that ~ 64% of the non-merger pairs are due to the projection 
effect, and ~ 36% of them are really 3D close pairs but fail to 
merge. The projection effect appears to be the main factor for 
close-pair contamination. It seems that it is likely high relative 
velocities to prevent 3D-close pairs from merging. 



sive host halos tend to be have higher values of (1 + 5 n ). 
However, total number of galaxies in massive clusters is 
smaller than that in groups or small clusters. In other 
words, high density (1 + 6 n ) is dominated by galaxies 
from groups or small clusters. 

(2) Although the pair fraction iV c among the SAMs 
displays large discrepancies among them and slightly de- 
viate from the observational data, their redshift evolution 
is consistently flat and in agreement with the same trend 
as the observati o nal re sults and as the theoretical results 
of IBerrier et"aTl (pOOl . 

(3) The normalized merger rate F mg and the volumet- 
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Fig. 11. — Top: The fraction of pairs with mass ratio fj, > 1/2 
(red), 1/3 (green), 1/4 (blue), and 1/6 (pink) to total pairs as a 
function of M% at two redshift range, 0.5 < z < 0.8 and 0.8 < z 
< 1.0. The fraction shows roughly no environmental dependence 
implying that the fi composition in pairs is roughly the same no 
matter what environment they reside. Bottom: The galaxy merger 
rate F mg as a function of Mj, for different mass ratio criterions at 
redshift 0.8 < z < 1.0. The similar shape among F mg profiles again 
infers that the impact due to mass ratio fi is not significant. 



slope is flatter in BowerOG and Font08 with a ~ 0.3 and 
steeper in DcLucia06 and Bertone07 with a ~ 0.6 — 0.7. 
That is, the merger rate between the low and high den- 
sities can differ by over an order of magnitude in cer- 
tain SAMs. By contrast, the profile of merger rate as 
a function of halo mass is not power laws but exhibits 
more complicated dependence. The group environment 
turns out to be the regime wher e galaxies merge most 
frequently, consistent with what iTran et ahl (|2008T ) ob- 
served; high merger rate in high local density regions is, 
in fact, due to significant contributions from galaxies in 
group environments. 

(6) The merger time T mg shows approximately no en- 
vironmental and no model dependence. The difference 
between the densest and under-dense regions is ~ 10%. 
Different environments have little influence on the merger 
timescale if these projected close pairs are true bound 
pairs. The merger time is found ~ 2 Gyr at z ~ 1. 

(7) Unlike T mg , merger probability of close pairs C mg 
depends on environments and on SAMs. Despite the dis- 
crepancies among different SAMs, their profiles all have 



ric merger rate R m g are both below the observational a similar trends that C mg drops as {\-\-8 n ) increase, and 



results. This is mainly because the ratio of the merger 
probability of close pairs C mg to the merger time scales 
T rn g (the normalization constant) found in the SAMs is 
much smaller than those adopted in observations. 

(4) When the pair fraction N c is expressed in terms of 
environment, N c is found to be higher in higher density 
or more massive Mh environment and the environmental 
dependence of N c evolves little with time. In addition, 
difference among different SAMs is found mainly in high 
density bins. 

(5) We observe strong environmental dependence of the 
galaxy merger rate in the SAMs, where higher density re- 
gions have larger merger rates. When the galaxy merger 
rate is expressed as F rng oc (1 + S n ) a , the logarithmic 



on the other hand, all SAMs have a peak in the bin with 
logio(Mh/M@) between 12 and 13. 

(8) Through evaluating the ratio of real 3D close pairs 
to 2D close pairs, we find that the projection effect is 
responsible for the environmental dependence of C mg . 
At z ~ 1, only ~ 31% of projected close pairs will 
merge eventually by 2 = 0. For those close pairs that 
do not merge, the projection effect appears to be the 
origin of major contamination. 

(9) Because of no dependence on environments for the 
fraction of pair satisfying the mass ratio criterion to total 
pairs and the similarity among different F mg profiles, we 
conclude that the impact of mass ratio on the merger 
rate in different environments should be a minor effect. 
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